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USER'S MANUAL FOR XTRAN2L (VERSION 1.2): A PROGRAM FOR SOLVING 

THE GENERAL-FREQUENCY UNSTEADY TRANSONIC 
SMALL-DISTURBANCE EQUATION 


David A. Seidel and John T. Batina 
NASA Langley Research Center 
Hampton, VA 23665 

Summary 

The development, use, and operation of the XTRAN2L program that solves 
the two-dimensional unsteady transonic smal 1 -di sturbance potential equation 
are described. The XTRAN2L program is used to calculate steady and unsteady 
transonic flow fields about airfoils and is capable of performing self- 
contained transonic flutter calculations. Operation of the XTRAN2L code is 
described, and tables defining all input variables, including default values, 
are presented. Sample cases that use various program options are shown to 
illustrate operation of XTRAN2L. Computer listings containing input and 
selected output are included as an aid to the user. 

Introduction 

Considerable research is presently being conducted to develop 
computational methods for calculating transonic unsteady aerodynamics for 
aeroelastic applications. The methods are being developed to provide accurate 
yet cost efficient means for calculating unsteady transonic airloads for the 
prediction of aeroelastic characteristics such as flutter and divergence. 

Computational developments and advances that are possible because of 
improvements made in computer technology and numerical solution techniques are 
described in Ref. 1. 

Solutions for unsteady transonic flows about oscillating airfoils have 
been made possible using computer codes such as LTRAN2.2 The original 
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LTRAN2 code was developed to compute time-accurate solutions of the low- 
frequency approximation of the transonic small -disturbance (TSD) equation. 
Houwink and Van der Vooren^ extended the range of applicability of the 
LTRAN2 code by incorporating the time-derivative terms in the airfoil and wake 
boundary conditions. The resulting code was termed LTRAN2-NLR. The XTRAN2L 
code, ^ developed at NASA Langley Research Center, is an extensive 
modification of LTRAN2-NLR that includes monotone differencing, nonreflecting 
farfield boundary conditions, an improved computational grid, a pulse 
transient capability, and an aeroelastic transient capability. Details of the 
algorithm used in XTRAN2L and the modifications made to LTRAN2-NLR are given 
by Whitlow. 4 Details of the grid development and pulse transient capability 
are given by Seidel, Bennett, and Whitlow. 5 The development of the 
aeroelastic transient capability was reported by Edwards, Bennett, Whitlow, 
and Seidel . 6 

The XTRAN2L program is operational at the NASA Langley Research Center on 
a CYBER-175 computer using FORTRAN IV. The program uses the OVERLAY 
capability to segment the program into 3 parts and reduce the field length 
required to load and run the program. The program requires 120 kbytes to load 
and run in overlayed form. On the CYBER-175, the program requires 0.25 CPU 
seconds per time step. 

The purpose of the present report is: (1) to summarize the XTRAN2L 

developments of Refs. 4-6, and (2) to describe the operation of the XTRAN2L 
(Version 1.2) code. The report reiterates the governing equations and program 
input descriptions given in the LTRAN2-NLR user's manual, ? and extends the 
information to include the XTRAN2L modifications and improvements. A brief 
description of XTRAN2L is presented first, highlighting the relevant equations 
and computational procedures. Operation of XTRAN2L is described next, along 
with tables defining all input variables, including default values. Finally, 
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application of XTRAN2L is illustrated by sample calculations for the NACA 
64A010A airfoil at a freestream Mach number of 0.78 and a mean angle of attack 
of 1.0°. In these calculations, various program options are used. Computer 
listings containing input and selected output are included as an aid to the 
user. 

Symbols 

a pulse amplitude in Eqs. (8a), (8b) and (8c) 

a^ distance in semi chords from midchord to elastic axis 

a n harmonic coefficient in Eq. (7) 

A coefficient defined in Eq. (2a) 

b semichord, c/2 

b n harmonic coefficient in Eqs. (6a) and (6b) 

B coefficient defined in Eq. (2b) 

c chordlength 

c a axial force coefficient 

c n normal force coefficient, positive up 

c n ^ control surface normal force coefficient, positive up 

c m pitching moment coefficient about pitch axis x r , positive nose 

down 

c m ^ control surface moment coefficient about hinge line x^, 
positive trailing-edge up 

c^j distance in semichords from midchord to hinge axis 

C coefficient defined in Eq. (2c) 

Cp pressure coefficient 

C Pu upper surface pressure coefficient 

C Pi , lower surface pressure coefficient 

★ 

Cp critical pressure coefficient 


i 
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f instantaneous position of airfoil, defined in Eqs. (4a) and (4b) 

f 0 (x) airfoil contour 

g(x,t) chord deformation 

g^(x) amplitude of harmonic chord deformation 

h(t) plunge displacement of elastic axis nondimensional ized by semichord, 
positive down (used for aeroelastic motion) 

H(t ) function for forced harmonic motion 

k reduced frequency, wb/U 

[k] stiffness matrix, defined in Eq. (10b) 

m airfoil mass per unit span 

M local Mach number 

Moo freestream Mach number 

[M] mass matrix, defined in Eq. (10a) 

p(x,t) arbitrary motion specified by the user 

r a radius of gyration of airfoil about elastic axis 

r^ radius of gyration of control surface about hinge line 

rj parameter to turn on/off time dependent term in airfoil boundary 

condition, defined in Eq. (3) 

r 2 parameter to turn on/off time dependent term in airfoil wake 

condition, defined in Eq. (12) 

rg parameter to turn on/off time dependent term in Cp and Mach 

number formulas, defined in Eqs. (14) and (16) 

t scaled time, kT 

£ nondimensional time for print and plot output, kx 

T time, (sec.) 

T c time at pulse center in Eqs. (8a), (8b) and (8c) 

U freestream speed 

Up flutter speed 

(u| aerodynamic load vector for aeroelastic motion, defined in Eq. (lOd) 
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w parameter that controls pulse width in Eqs. (8a), (8b) and (8c) 

x streamwise coordinate relative to leading edge, measured in 

chordlength 

xf control surface leading edge location 

xh control surface hinge axis location 

x r pitch axis location 

x a distance in semichords from elastic axis to center of mass 

distance in semichords from hinge axis to control surface center 
of mass 

{x} displacement vector for aeroelastic motion, defined in Eq. (10c) 

z coordinate normal to freestream, measured in chordlength, 

positive up 

Z(t) nondimensional plunge displacement, positive up (used for forced 
harmonic and pulse motions) 

Zj harmonic plunge amplitude 

a(t ) pitch motion, positive nose up 

a 0 mean angle of attack 

oq harmonic pitch amplitude 

p(t) control surface motion, positive trailing edge down 
3 0 control surface mean deflection angle 

control surface harmonic amplitude 
y ratio of specific heats 

ACp lifting pressure coefficient 

aT time step size 

6 airfoil maximum thickness-to-chord ratio 

[ej integral of state-transition matrix, used for aeroelastic motion 

m airfoil mass ratio, m/irpab? 

p density 

Pa. freestream density 
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t nondimensional time, UT/b 

<f> disturbance velocity potential 

[<t] state-transition matrix, used for aeroelastic motion 
a) angular frequency 

0 )^ uncoupled natural frequency of plunge, rad/sec 

oo a uncoupled natural frequency of pitch, rad/sec 

uncoupled natural frequency of control surface deflection, rad/sec 

Description of XTRAN2L 

The XTRAN2L code is used to calculate time-accurate, finite-difference 
solutions of the nonlinear, small -disturbance, potential equation for 
two-dimensional transonic flow. The code can be used to calculate steady and 
unsteady transonic flow fields about airfoils and is capable of treating 
forced harmonic, pulse, or aeroelastic transient type motions. In this 
section, the XTRAN2L code is described including the relevant equations and 
computational procedures. 


Transonic Small -Disturbance Equation 

The flow field is described by the complete unsteady transonic small - 
disturbance equation 


A *tt + B Vt = C *xx + * 


zz 


where $ is a disturbance velocity potential scaled by c[6f*r (y* + 1)] 


( 1 ) 

2/3 


2 

Y* = 2 - (2 - and the spatial coordinates x and z and time t are 

P 1 / O 1 

normalized from the physical quantities by c, c/[6M^ (y* + 1)] ' and k , 
respectively. The coefficients A, B, and C are defined as 
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(2a) 


A = c 2 k 2 Mf/U 2 [6M 2 (y* + 1)] 2/3 
8 = 2 c k M*/U [6M^ (y* + 1)] 2/3 (2b) 

C = (1 - M 2 )/[6M^ (y* + 1)] 2/3 - ( y* + 1)4> X (2c) 

A time-accurate alternating-direction implicit (ADI) finite-difference scheme 
is used to solve Eq. (1). Details of the ADI algorithm are given by 
Whitlow. A 


Airfoil Boundary Condition 

The airfoil flow tangency boundary condi tion? is 



+ r, 


f 


t 


( 3 ) 


where 

f = f 0 (x) + g(x,t) + Z(t) - a(t) [x - x p ] + p(x,t) ; x < x f (4a) 


and 

f = f 0 (x) + g(x,t) + Z(t) - a(t) [x - x p ] + p(x,t) 

- B(t) [x - x h ] ; x > x f (4b) 

In Eq. (3), rj is a parameter specified by the user to be either 0 or 1. 
The low frequency airfoil boundary condition is selected by setting rj to be 
equal to zero. In Eqs. (4a) and (4b), f 0 (x) describes the airfoil contour, 
g(x,t) is the chord deformation, and Z(t), a(t), and B(t) are the 

nondimensional plunge displacement, pitch, and trailing edge control surface 
deflection, respectively. Three types of motion may be analyzed using XTRAN2L 
as described below. The motion types are forced harmonic, pulse, and 
aeroelastic. The capability for specifying an arbitrary motion is provided 
through the use of the function p(x,t). 
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Forced Harmonic Motion. - Typically, unsteady aerodynamic forces are 
determined by calculating several cycles of forced harmonic oscillation with 
the last cycle providing the estimate of the forces. Following Ref. 7, the 
equations describing forced harmonic motion are 


g(x,t) = g^x) H(t) (5a) 

Z(t) = Z x H(t) (5b) 

a(t) = a Q + ot 1 H(t) (5c) 

»(t) = + 3 1 H(t) (5d) 


where a 0 and tJ 0 are the airfoil mean angle of attack and control surface 
mean deflection angle, respectively; gi(x) is the amplitude of harmonic 
chord deformation, and Zj, c^, and 01 are the amplitudes of harmonic 
plunge, pitch, and control surface motions, respectively. Two alternatives 
exist for the definition of the harmonic function H(t). The first definition 
is 

N 9 

H(t) = l b sin^ nkx ; 0 < t < tt (6a) 

n=l n 2nk 

N 

H(t) = I b n sin nkx ; x > (6b) 

n=l 

where N is the number of harmonics to be selected by the user. Equation (6a) 
is used to provide a smooth start over the first quarter period of each 
harmonic which reduces the starting transient. The second definition for H(t) 
is 

N 

H(t) = l a cos nkx ; x>^0 (7) 

n=l n 
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Real and imaginary parts of the unsteady aerodynamic forces in the 
frequency domain may be determined two ways. First, the method used in the 
LTRAN2-NLR code was retained. This method calculates the coefficients of the 
first harmonics of the airloads based on information from each quarter period 
of the last cycle of harmonic motion. The procedure assumes that the time 
histories of the normal force and moment coefficients are periodic and simple 
harmonic. The relevant equations are presented in Ref. 7. In the second 

procedure, Fourier analyses of the last cycle of motion are performed to 
determine all relevant harmonic components. This procedure is the preferred 
method since the only assumption is that of periodicity. 


Pulse Motion . - As an alternative to calculating unsteady aerodynamic 
forces by forced harmonic oscillation, the unsteady forces may be obtained 
indirectly from the response due to a smoothly varying, exponentially shaped 
pulse. 5 In this procedure, the airfoil is given a small prescribed pulse in 
a given mode of motion and the aerodynamic transients are calculated. For 
plunge motion, the input pulse is given by 

Z(t) = ae" w " T / (8a) 

for pitch motion, the input pulse is given by 


a(t) 


+ ae 


-w (T - 




(8b) 


and for control surface motion, the input pulse is given by 


ii(t) = tJ + ae" w ^ " V 


(8c) 


where a is the pulse amplitude, w is a parameter that controls the pulse 
width, and T c is the time at the pulse center. The pulse width is inversely 
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proportional to w, such that as the value of w increases, the width of the 
pulse decreases. To ensure adequate resolution of harmonic components up to a 
reduced frequency of k = 2.0, the product of the chosen values of w and aT^ 
should be near the product of the default values. The harmonic response is 
obtained by a transfer-function analysis using fast Fourier transforms (FFT). 
Specifically, the aerodynamic forces in the frequency domain are determined by 
dividing the FFT of the time histories of the output normal force and moment 
coefficients by the FFT of the input airfoil pulse motion. The calculation is 
performed external to the XTRAN2L code as a separate job step and thus the 
post-processing program is not included within XTRAN2L. Use of the pulse 
transfer-function analysis gives considerable detail in the frequency domain 
with a significant reduction in cost over the alternative method of 
calculating multiple oscillatory responses, as shown in Fig. 13 of Ref. 5. As 
an example. Fig. 1 shows a pitch pulse input (Fig. 1(a)), a typical normal 
force coefficient output (Fig. 1(b)), and the resulting unsteady forces (Fig. 

1(c)). 


Aeroelastic Motion . - An aeroelastic transient capability is available 
within XTRAN2L wherein the structural equations of motion are coupled with the 
aerodynamic solution procedure for simultaneous time-integration. The 
equations of motion for a three degree-of-f reedom (plunge, pitch, and control 
surface motions) aeroelastic system may be written as 


[M] ix} * LK] tx) = ij $ |u| 


( 9 ) 


where [m], the mass matrix, [k], the stiffness matrix, {x}, the displacement 
vector, and {u } , the aerodynamic load vector, are defined as 
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(10a) 


(10b) 




The dot (•) above a variable denotes differentiation with respect to time T, 
and the structural parameters appearing in Eqs. (9) and (10) are defined in 
Fig. 2. In the aerodynamic load vector, the steady-state coefficients 
c n 0 » Cm 0 , and c^ are subtracted from the total coefficients. 
Thus the aeroelastic response (i.e., h(t), a(t), and 3(t)) represents motion 


11 



about the steady-state condition. Aeroelastic equations of motion for any of 
the three possible two degree-of-freedom systems are determined by deleting 
the appropriate rows and columns in Eqs. (9) and (10). 

The aeroelastic response is determined by numerically integrating the 
equations of motion in time using a structural integrator based on the 
state-transition matrix. 6 More specifically, the time-integration of Eq. 
(9) is performed by advancing the aeroelastic solution from time level n to 
time level n+1 using an algorithm described by 



[#] 






T r- 


J 


M 


n 



(ID 


where [<&J and [aj are the state-transition matrix and its integral matrix, 
respectively. These matrices are functions of the aeroelastic parameters and 
are computed within the XTRAN2L code. Details of the aeroelastic solution 
procedure may be found in Ref. 6. In general, several aeroelastic responses 
are calculated for a range of freestream speeds, U, at a given Mach number. 
The values for U are selected such that stable (converging) and unstable 
(diverging) responses result. The flutter speed Up, where the damping is 
zero, is calculated by interpolating the critical damping of the responses. 
Damping of the aeroelastic transients is typically determined using the modal 
identification program of Bennett and Desmarais.8 The program uses complex 
exponential functions and a least-squares curve-fitting technique to estimate 
damping, frequency, magnitude, and phase of the aeroelastic modes from the 
transient responses. This calculation is performed as a separate 

post-processing job step and thus the program of Ref. 8 is not included within 
the XTRAN2L code. 


12 



Wake Condition 


The wake is represented as a slit downstream of the airfoil trailing 
edge. Along the slit the condition? 

[\] + r 2 [* t ] = 0 ; z = 0, x > 1 (12) 

is imposed, where [ ] denotes the jump in the indicated quantity across the 
wake. The parameter r 2 is specified by the user to be either 0 or 1. The 
low frequency version of the wake condition is specified by setting r 2 to be 
equal to zero. 

Nonreflecting Farfield Boundary Conditions 

The LTRAN2 and LTRAN2-NLR codes use steady -state farfield boundary 
conditions which cause disturbances to be reflected from the computational 
boundaries. The boundaries were subsequently located far away to minimize the 
possibility of disturbances reflecting back into the flow field and 
contaminating the solution. Therefore, nonreflecting farfield boundary 
conditions were developed and implemented in XTRAN2L by Whitlow. 4 These 
boundary conditions, consistent with the TSD equation (Eq. (1)), are 
summarized in Fig. 3, taken from Ref. 4. 


Pressure Coefficient and Mach Number 

The pressure coefficient C p , may be calculated using either of two 
formulae. 7 The exact expression is given by 



where 


1 

p = [1 - Ijjr M l ( 2<J> X + 4> x + *2 + Y_1 


(13b) 
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The linearized TSD approximation is given by 


C p " - 2 K 


+ r lW 


( 14 ) 


The local Mach number M, may also be calculated by either of two 
formulae. 7 The exact expression is 


M 


2 



d + *x>‘ 


.Y-l 



(15) 


and the linearized TSD approximation is 

M 2 = M 2 [1 + (y* + 1)4> X + r 1 r 3 ( Y - 1)M f <t> t ] (16) 

The user may also select low-frequency versions of the linearized 
equations, Eq. (14) and Eq. (16), whereby the time-derivative terms are 
omitted. This is done by defining the product r^ to be equal to zero. 
Otherwise the product r^ is set equal to unity. 


Normal Force, Axial Force, and Pitching Moment Coefficients 

For the airfoil, the normal force coefficient c n , axial force 
coefficient c a , and pitching moment coefficient c m , are defined as 


fl 


c n = 


Jo 


AC dx 
P 


(17a) 


c 


a 





f ) dx 
x „ ' 


(17b) 


'1 

c = - AC (x - x J dx 
m p ' r' 

JO 


(17c) 
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where ACp = Cp u - Cp^, and the subscripts u and % denote the upper 
and lower airfoil surfaces, respectively. 

For the control surface, the normal force coefficient c nf , and 
pitching moment coefficient c m ^ are defined as 



(\ dx 


(18a) 



(" AC p (x - x h ) dx 


(18b) 


By substituting Eq. (14) for the pressure coefficients appearing in Eqs. 
(17) and (18), and evaluating several of the integrals by parts leads to 


c n ' 2 4 *TE + 2r l r 3 f 4 *t dx 

Jo 


(19a) 


c a = iV \ ~ C D f x ) dx 

a 1 p u x u P £ x £ 


(19b) 


0 


fl 


c m ■ 2 0 - x r> 4 *TE + 2 


A<}> dx + 2r^r 3 


A4> t (x - x p ) dx 


(19c) 


for the airfoil, and 


c n f ■ 2 4 »TE - 2 4 »x f + 2r l r 3 f 4 *t dx 


(20a) 


c - 2 (1 - x h ) A4 > je - 2 (x f - x h ) A* + 2 | A<j> dx 


•1 

+ 2 r 1 r 3 | A4> t (x - x h ) dx 
x r 


(20b) 
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for the control surface. In Eqs. (19) and (20), a<|> = $ u - <p a , and a<j»je 

and A<t> Xir are the values of a<i> at the airfoil trailing edge and control 
surface leading edge, respectively. In XTRAN2L, the normal force, axial 
force, and pitching moment coefficients are calculated using Eqs. (19) and 
(20) when the linearized equation is selected for the pressure coefficient 
(Eq. (14)). This is computationally more accurate than integrating the 
lifting pressure coefficients as indicated in Eqs. (17) and (18). Note that 
if a particular case has no control surface then Xf = xh = 1 leads to 
c nf = 0 and c mf . = 0. 


Operation of XTRAN2L 
General Information and Input Data 

In this section, general information is given on the operation of XTRAN2L 
including input data and recommended operating procedures. Data to define 
the problem and control program execution are defined using six namelists. 
The input data consists of a single TITLE card followed by the six namelists 
in the following order: 


AEROEL - Definition 
INPUT - Definition 
MESH - Definition 
AMPLI - Definition 
ORD - Definition 
COND - Definition 


of airfoil motion and aeroelastic parameters 
of program solution parameters and constants 
of user-input extended fine grid 
of periodic time-dependent airfoil motion 
of airfoil geometry 
of airfoil chord deformation 


Namelists INPUT, MESH, AMPLI, ORD, and COND have been retained from the 
LTRAN2-NLR code. A new namelist, AEROEL, has been added to allow for 
specification of the type of airfoil motion (forced harmonic, pulse. 
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aeroelastic, or no motion) and for input of pulse or aeroelastic parameter 
values. Tables 1-6 list the parameters input in each namelist and give a 
short description of each parameter including its default value. The tables 
are presented in the same format as that of the LTRAN2-NLR user's manual. 7 
Any parameter not explicitly defined in a NAMELIST statement assumes its 
default value. 

To start the solution of the unsteady flow equation, the program requires 
an initial, steady-state solution for the flow field at T = 0. The 
steady-state solution can be calculated in one of three ways which are 
selected using the IREAD parameter in namelist INPUT. The first method is to 
define the initial flow field as an undisturbed stream (IREAD = 0). The 
second method is to solve for a steady-state solution using a successive 
line-overrelaxation (SLOR) scheme (IREAD = 1). The third method is to solve 
for a steady-state solution using the unsteady TSD scheme, termed the 
"unsteady solver", with no airfoil motion (IREAD = 0 and IRESP = 5). Once a 
converged steady-state solution has been achieved, the program is then 
restarted using the converged solution as the initial condition (IREAD = 5). 
The SLOR solution procedure normally involves using a sequence of three 
default grids defined as coarse, fine, and extended fine, to accelerate 
convergence. The default extended fine grid used in XTRAN2L was developed to 
increase the accuracy of the unsteady solution while decreasing the number of 
grid points used in the LTRAN2-NLR code. 5 Grid point locations for the 
default extended fine grid are listed in Table 7. The grid distribution in 
the vicinity of the airfoil is shown in Fig. 4. The coarse and fine grids 
were defined to be subsets of the extended fine grid. The user may choose to 
use the default extended fine grid or input a new extended fine grid using 
namelist MESH. 
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Calculating steady solutions using the unsteady solver requires starting 
from an initial undisturbed stream and solving the unsteady TSD equation for 
no airfoil motion until a converged steady-state solution is reached. All 

computations using the unsteady solver are performed on the extended fine grid 
only. Experience has shown that for most cases the unsteady solver converges 
to a steady-state solution as fast or faster than the SLOR algorithm. Use of 
the unsteady solver is the preferred approach since for lifting cases the two 
algorithms converge to slightly different steady-state solutions. 5 For 

consistency with the unsteady solution it is therefore recommended that the 
steady -state starting solution be calculated using the unsteady solver. 

The unsteady solver can solve three combinations of the TSD equation and 
boundary conditions. These are denoted as the LTRAN, HYTRAN, and EXTRAN 
options. The equation set to be solved is chosen using the IEQFLAG parameter 
in namelist AEROEL. Selecting the LTRAN option gives the low-frequency 
approximation to the complete unsteady TSD equation (Eq. (1)) by defining the 
coefficient "A" in the equation to be zero. In addition, the time-derivative 
terms in the airfoil and wake boundary conditions (Eqs. (3) and (12), 
respectively), and in the pressure and local Mach number equations (Eqs. (14) 
and (16), respectively), are automatically deleted by the program setting the 
parameters rj and (variables IR and IMODFR in namelist INPUT) to be 
equal to zero. Selecting the HYTRAN option defines the coefficient "A" in 
Eq. (1) to be zero. The time derivative terms are included in the airfoil and 

wake boundary condition equations and in the pressure and local Mach number 

equations by setting r^ and r 2 to be equal to unity. This option is the 
capability implemented in LTRAN2-NLR.? In choosing the LTRAN and HYTRAN 
options, the farfield boundary conditions in Fig. 3 reduce to a set consistent 
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with the low frequency approximation. Selecting the EXTRAN option defines the 
coefficient "A" in Eq. (1) to be as described in Eq. (2a). The time 
derivative terms are also included in the associated boundary conditions. 

The unsteady solver calculates the time varying flow field by employing a 
time marching solution procedure. For oscillatory cases, the time step is 
implicitly defined by the user based upon the reduced frequency k, scaled 
velocity U/c, and number of time steps per cycle NSPC, as 

, T if 

NSPC*k*U/c (21) 

For cases other than oscillatory motion, the time step is explicitly defined 
by the user. In these cases, the reduced frequency k is implicitly defined 

based on AT, U/c, and NSPC, using Eq. (21). 

For both steady and unsteady cases, the airfoil contour is defined in 

terms of the upper and lower surface slopes at the grid points, input in 

namelist ORD. The program will interpolate airfoil slopes which are not 

defined at the grid points and thus calculate new slopes, but this approach is 
not recommended. Before running the XTRAN2L program, a pre-processor is 
normally used to convert the airfoil geometry into slopes and displacements at 
the grid points. This allows the user to view the geometry to be input to the 
program and to take any required corrective action. This is especially true 
for measured airfoil ordinates where errors or inaccuracies in the ordinates 
may cause unrealistic airfoil slopes to be calculated. 

The program has a restart capability designed for continuing the steady 
SLOR solution or the unsteady solver from a previously saved flow field 
solution. The restart capability is typically used to restart the SLOR 
solution to achieve better steady-state convergence (IREAD = 3 or 4), to start 
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an unsteady computation from a previously calculated steady-state solution 
(IREAD = 5), or to continue an unsteady computation further in time 

(IREAD = 5). 

Files and Logical Units 

The XTRAN2L program uses the logical units TAPE3, TAPE4, TAPE5, TAPE6, 
TAPE7, and TAPE8. The logical units TAPE5 and TAPE6 are used for program 
input and output, respectively. The logical units TAPE3 and TAPE4 are used 
for restart files. The data to be saved for a restart is written to TAPE4 and 
a restart reads the required data from TAPE3. The data is in unformatted 
binary form and is input and output using BUFFER IN and BUFFER OUT 
statements. The BUFFER IN and BUFFER OUT statements are called in subroutines 
GETS and PUTS, respectively. The call to GETS and PUTS is of the form: 

CALL GETS (INDEX, ARRAY, NELEMNT) 


where , 

INDEX - logical unit number 

ARRAY - array of data to be input/output 

NELEMNT - number of data elements to be input/output 

The data is written to TAPE4 using the following FORTRAN statements: 

CALL PUTS (4,BUF1,2) 

CALL PUTS (4,BUF2,40) 

CALL PUTS (4,BUF3,5) 

CALL PUTS (4,VU0,80) 

CALL PUTS (4, VL0,80) 
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CALL PUTS (4,PTU(1,ITM1), JMAX1) 

CALL PUTS (4,P0U,JMAX1) 

CALL PUTS ( 4 , PTL ( 1 , ITM1 ) , JMAX1 ) 

CALL PUTS (4,P0L,JMAX1) 

DO 704 J=l, JMAX1 
704 CALL PUTS (4,P0(1,J),LMAX1) 

CALL PUTS (4,PTT(1,1),4880) 

CALL PUTS ( 4 , PTTU ( 1 ) , 80 ) 

CALL PUTS (4,PTTL(1),80) 

where, 

BUF1 - Array containing circulation value and time from last time step 

BUF2 - Array containing array sizes, airfoil aeroelastic response 

from last time step and calculated forces from last time step 

BUF3 - Array containing calculated forces from last time step 

VUO/VLO - Array of airfoil upper/lower surface boundary condition 
values from last time step 

PTU/PTL - Array of airfoil upper/lower surface potentials from the next 
to last time step 

JMAX1/LMAX1 - Number of points in extended fine grid in x/z, respectively 

POU/POL - Array of airfoil upper/lower surface potentials from last 
time step 

PO - Array of flow field potentials from last time step 

PTT - Array of flow field potentials from next to last time step 

PTTU/PTTL - Array of airfoil upper/lower surface potentials from next to 
last time step 

The logical units TAPE7 and TAPE8 are used for storing results at each time 
step for post-processing. If the IBUF parameter in namelist INPUT is set 
equal to 1, the program outputs time histories of aeroelastic quantities, 
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pressure coefficients, and forces in unformatted binary to TAPE8. For 
restarts, TAPE7 is provided to allow the merging of previous time histories 
from a TAPE8 file onto a new file TAPE8. The program copies TAPE? onto TAPE8 
and then proceeds as normal. This option is not recommended due to the large 
size required to store a TAPE8 file. A blank TAPE7 file causes the program to 
skip the copying step on restarts and continues with normal execution. The 
data is written to TAPES using the following FORTRAN statements: 

CALL PUTS (8, TITLE, 8) 

CALL PUTS (8, 1 INPUT, 49) 

CALL PUTS (8,X,NPTSARF) 

CALL PUTS (8,FXU ,NPTSARF ) 

CALL PUTS (8,FXL,NPTSARF) 

00 10 IA=1, NITER 
CALL PUTS (8,BUF4,9) 

CALL PUTS (8,CPU,NPTSARF ) 

CALL PUTS (8,CPL,NPTSARF ) 

CALL PUTS (8,XN,6) 

10 CONTINUE 

where, 

TITLE - array containing 80 character title 

IINPUT - array containing program variables and constants 

X - array containing x coordinates of grid points on airfoil 

/ 

NPTSARF - number of x-grid points on airfoil 

FXU/FXL - array containing airfoil upper/lower surface 
displacements at x-grid points 
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NITER - number of iterations of time history data written 
=1 - steady solver solution only 
=MAXIT - unsteady solution only, no restart 
=MAXIT+1 - unsteady solution started from steady solver or 
restart file 

BUF4 - array containing nondimensional time and calculated force 
coefficients 

CPU/CPL - array containing upper/lower surface pressure coefficients 
on airfoil 

XN - array containing aeroelastic displacements of airfoil 


The do-loop illustrated above is executed at every time step. If a restart is 
performed, the first set of data written by the do-loop is for the data read 
from the restart file. 


Sample Calculations 

Sample calculations were performed for the NACA 64A010A airfoil 9 (NASA 
Ames model) to illustrate various program options. Four cases were run to 
demonstrate the program's capabilities. The first case is a steady-state 
solution using the unsteady solver. The next three cases are a simple 
harmonic pitching motion of the airfoil, an airfoil pitch pulse which is used 
to calculate the generalized aerodynamic forces, and a transient aeroelastic 
response of the airfoil at a condition just above the flutter speed. The last 
three cases are restarted from the steady-state solution. 

Steady Calculations 

Results from the steady calculation are listed in Appendix A. The input 
file read by the program is listed first, followed by the program output. 
For brevity, the output included in Appendix A is only a partial listing of 
the output generated by the program for the case run. The output shown 
includes the banner page, echo of input parameters, airfoil description. 
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computational grid listing and output of pressure distributions and forces 
calculated at selected iterations. The unsteady solver was used to converge 
to a steady flow field solution at M«, = 0.78 and oq = 1.0° starting from 
freestream initial conditions with no airfoil motion. The calculation was run 
for 1024 iterations to ensure convergence to a steady-state solution. Airfoil 
pressures and forces are printed every 128 iterations. The pressure dis- 
tributions and forces listed in Appendix A are for the first and last 
iterations for which calculations were performed. The time listed with each 

set of pressure output is a nondimensional time based on harmonic motion, E. 

For harmonic motion, l is the measure of the angular position in radians in 
the harmonic cycle. Also output is the equivalent angular position in 
degrees. For other than harmonic motion, the time output can be considered as 
a nondimensional time. The time histories of airfoil force coefficients and 
the final pressure distributions are shown in Fig. 5. The plots were created 
by a post-processing program which reads the TAPE8 file and uses the 
information to generate appropriate plots for the case run. 

Unsteady Calculations 

Forced Harmonic Motion . - Results from a forced harmonic oscillation of 
the airfoil in pitch are listed in Appendix Bl. The input file is listed 
first, followed by selected program output. The airfoil was oscillated for 
four cycles about the quarter-chord at a reduced frequency of k = 0.075 with 
an amplitude of = 0.5° about the mean angle of attack oq = 1.0°. The 
flow field was calculated in 1° increments of a cycle of motion. Included in 
the output are the pressure distributions for the first and last iterations 
and time histories of the airfoil forces for the last cycle of oscillation. 
In the output, the pitching moment coefficient, normal force coefficient, and 
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control surface moment coefficient are labeled as "CM", "CN", and "CMF", 

respectively. A Fourier analysis of the forces is performed using the last 
cycle of airfoil motion and the first nine harmonics are output. The real and 
imaginary components of the forces, labeled as "RE" and "IM", respectively, 
are scaled by the amplitude of the motion (in radians for both airfoil pitch 
and control surface motions). Also included in the harmonic quantities output 
is a measure of the amount of higher harmonic components in the force time 
histories, labeled as "ECM", "ECN", and "ECMF". The quantities are defined as 
the magnitude of the nth harmonic component divided by the magnitude of the 
fundamental, or first harmonic component. The zeroth harmonic is defined as 
the mean or average value. 

The first five harmonics of the airfoil pressures calculated for the last 
cycle of motion are included in Appendix B1 for reference. The pressure 
harmonics are calculated by a post-processing program using the information 
stored on TAPE8. The time histories of the airfoil motion and corresponding 
force coefficients are plotted in Figs. 6(a) and 6(b), respectively. Figures 
6(c) -6 (f ) show the mean pressure distribution and the first harmonic pressure 
distributions calculated by the post-processor. 

Pulse Motion . - Results from an airfoil pitch pulse are listed in 
Appendix B2. The input file is listed first followed by selected program 
output. The airfoil was pitched about the quarter-chord with maximum 

amplitude of a = 0.5° from the mean angle of attack of u 0 = 1.0°. The 
exponential pulse was centered at iteration 30, and the unsteady solver was 
run until the flow field returned to the steady-state condition. It is 
important to ensure that the flow field returns to steady-state conditions in 
order to allow accurate calculation of the generalized aerodynamic forces. 
The pressure distributions for the first and last iteration are listed. 
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Included in the output are the generalized aerodynamic forces calculated by a 
post-processing program from the information contained on TAPE8. The time 
histories of the airfoil motion and resulting force coefficients are plotted 
in Figs. 7(a) and 7(b), respectively. Figures 7(c) and 7(d) show the 
generalized aerodynamic forces calculated by the post-processor. 

Aeroelastic Motion . - Results from aeroelastic motion of the airfoil are 
listed in Appendix B3. The input file is listed first followed by selected 
program output. The aeroelastic parameters used were the program default 
values: ah = -2.0, x a = 1.8, r a = 1.865, y = 60.0, = 100.0 

rad/sec, and co a = 100.0 rad/sec. The default values are those of Case A of 
IsogailO which were chosen to have normal modes similar to those of a 
streamwise section near the tip of a sweptback wing. The wind-off coupled 
plunge and pitch frequencies are 71.3 and 535.7 rad/sec, respectively. The 
pivotal point for the plunge mode is located 1.44 chordlengths forward of the 
leading edge; the pivotal point for the pitch mode is located 0.068 
chordlengths ahead of midchord. For the sample aeroelastic calculation 

presented here, the scaled velocity U/c was set equal to 330.0, which is 
slightly above the flutter speed of the NACA 64A010A airfoil at M = 0.78 and 
«o = 1.0°. The program automatically calculates the state-transition matrix 
[$] and its integral matrix [ 0 ] in Eq. (11) from the aeroelastic quantities 
input in namelist AER0EL. The output from this calculation is included in 
Appendix B3. The [<&] matrix is labeled in the output as "PHI" and the 
"THETA*B-PRIME" matrix is defined as 

THETA*B-PRIME = ^ (^) 2 [ 0 ] jj°j-i| (22) 

The "I.C. VECTOR" listed in the output contains the initial values used for 
the aeroelastic displacement vector {x}. The pressure distributions output 
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for the first and last iterations are included in the Appendix. The time 
histories of the airfoil motion and corresponding force coefficients are 
plotted in Fig. 8. As shown, the airfoil motion is slightly divergent for the 
case selected. 


Acknowledgment 

The authors wish to acknowledge the National Aerospace Laboratory, NLR, 
the Netherlands, for the thorough presentation of equations and detailed input 
description for the LTRAN2-NLR code reported in their user's manual, 7 the 
organization and format of which were adopted in the present report. 

References 

iBallhaus, W. F., "Computational Aerodynamics and Supercomputers," NASA 
TM 85887, January 1984. 

^Ballhaus, W. F., and Goorjian, P. M., "Implicit Finite-Difference 
Computations of Unsteady Transonic Flows About Airfoils," AIAA Journal , 
Vol . 15, No. 12, December 1977, pp. 1728-1735. 

3 Houwink, R., and Van der Vooren, J., "Improved Version of LTRAN2 for 
Unsteady Transonic Flow Computations," AIAA Journal , Vol. 18, No. 8, 
August 1980, pp. 1008-1010. 

“^Whitlow, W., Jr., "XTRAN2L: A Program for Solving the General 

Frequency Unsteady Transonic Small Disturbance Equation," NASA TM 85723, 
November 1983. 

5 Seidel, D. A., Bennett, R. M. , and Whitlow, W. , Jr., "An Exploratory 
Study of Finite-Difference Grids for Transonic Unsteady Aerodynamics," AIAA 
Paper No. 83-0503, presented at the AIAA 21st Aerospace Sciences Meeting and 
Technical Display, Reno, NV, January 10-13, 1983, also NASA TM 84583, December 
1982. 


27 



^Edwards, J. W. , Bennett, R. M. , Whitlow, W., Jr., and Seidel, D. A., 
"Time-Marching Transonic Flutter Solutions Including Angle-of -Attack Effects," 
Journal of Aircraft , Vol. 12, No. 11, November 1983, pp. 899-906. 

?Anon., "Users Manual for LTRAN2-NLR; A Programme for the Calculation 
of Inviscid Transonic Flow About Thin Airfoils in Moderately Slow Unsteady 
Motion." NLR Memorandum W-79-003, June 1979. 

8 Bennett, R. M. , and Desmarais, R. N., "Curve Fitting of Aeroelastic 
Transient Response Data with Exponential Functions in Flutter Testing 
Techniques," NASA SP-415, 1975, pp. 43-58. 

^Bland, S. R., "AGARD Two-Dimensional Aeroelastic Configurations." 
AGARD AR-156 , August 1979. 

l°Isogai , K., "On the Transonic-Dip Mechanism of Flutter of a Sweptback 
Wing." AIAA Journal , Vol. 19, No. 7, July 1979, pp. 793-795. 


28 



Table 1. Description of input parameters to NAMELIST AEROEL. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

IRESP 


0 

Selects type of airfoil motion 

0 - Harmonic oscillation 

1 - Aeroelastic transient response 

2 - Plunge pulse transient response 

3 - Pitch pulse transient response 

4 - Control surface pulse transient 

response 

5 - No airfoil motion allowed 

NELAST 


1 

Selects type of aeroelastic motion 

1 - 2-DOF, plunge and pitch 

2 - 2-DOF, plunge and control surface 

3 - 2-DOF, pitch and control surface 

4 - 3-DOF, plunge pitch and control 

surface 

UC 

U/c 

100.0 

Scaled freestream velocity, used to 
define time step &T (IRESP = 0) or 
reduced frequency k (IRESP * 0), Eq . ( 21 ) 

XMU 

y 

60.0 

Airfoil mass ratio 

DELT 

AT 

.002618 

Time step (sec), calculated by program 
using Eq. (21) if IRESP = 0 

IEQFLAG 


3 

Selects equation set 

1 - LTRAN option 

2 - HYTRAN option 

3 - EXTRAN option 

AMPLTD 

a 

.017435 

Amplitude of pulse transient 

WDTHPLS 

w 

10000.0 

Exponential constant for pulse 

TZEROST 

T c /AT 

17.5 

Center of pulse (iterations) 

XALPHA 

X a 

1.8 

Nondimensional distance from elastic 
axis to airfoil center of mass 

XBETA 

x y 

.020833 

Nondimensional distance from hinge 
line to control surface center of mass 

RALPHA 

r 

a 

1.865 

Radius of gyration of airfoil about 
elastic axis 

RBETA 


.027169 

Radius of gyration of control surface 
about hinge line 
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Table 1. Concluded. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

OMEGAH 

w h 

100. 

Uncoupled natural frequency of plunge 

OMEGAA 

0) 

a 

100 . 

Uncoupled natural frequency of pitch 
about elastic axis 

OMEGAB 

"3 

500. 

Uncoupled natural frequency of control 
surface deflection about hinge line 

AH 

3 h 

-2. 

Nondimensional elastic axis location 

CBETA 

C 3 

.5 

Nondimensional control surface hinge 
location 

XIC(l) 

h(0) 

.01 

Initial plunge displacement 

XIC ( 2 ) 

a(0) 

0 . 

Initial pitch rotation 

XI C ( 3 ) 

3(0) 

0 . 

Initial control surface deflection 

XIC(4) 

h(0) 

0 . 

Initial plunge velocity 

XIC (5) 

a(0) 

0 . 

Initial pitch velocity 

XI C ( 6 ) 

3(0) 

0 . 

Initial control surface velocity 


3 


a 











Table 2. Description of input parameters to NAMELIST INPUT. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

XK 

2k 

.2 

Reduced frequency (based on chord), cal- 
culated by program using Eq. (21) if 
IRESP * 0 

XM 

M 

00 

.85 

Freestream Mach number 

GAM 

Y 

1.4 

Ratio of specific heats 

IREAD 


1 

Starting conditions 

0 - Unsteady computation starting from 

zero initial conditions 

1 - Compute steady-state initial 

conditions starting from zero, then 
do the unsteady computation 

2 - Compute steady-state solution 

starting from zero and stop 

3 - Read starting solution, compute 

steady-state solution and stop 

4 - Read starting solution and compute 

steady-state initial conditions, 
then do the unsteady computation 

5 - Read initial conditions and do the 

unsteady computation 
Options 3, 4, and 5 require a file 
which has been generated in a previous 
run 

ICNVRG 


1 

Selects method for converging steady- 
state solution 

0 - Converge solution on extended fine 

grid only 

1 - Converge to approximations to 

steady-state on coarse and fine 
spaced grids, then final solution on 
extended fine grid 

MAX IT 


360 

Number of time steps 

NSPC 


120 

Number of steps taken per cycle of 
oscillation (IRESP * 0), used to define 
time step AT (IRESP = 0) or reduced 
frequency k (IRESP * 0), Eq. (21) 
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Table 2. Continued. 


PARAMETER VARIABLE DEFAULT DESCRIPTION 

ITRAN 0 Selects airfoil motion 

0 - Pitch, plunge and control surface 

deflection 

1 - Pitch, plunge, control surface 

deflection and chord deformation 

2 - Arbitrary user-specified time- 

dependent motion only, user 
specifies p(x,t), Eqs. (4) 

ISICO 0 Selects method for calculating the time- 

dependent angle of attack, Eqs. (6) & (7) 

0 - Sinusoidal variation 

1 - Cosinusoidal variation 

ICPRNT 60 Time step print increment for inter- 

mediate Cp output 

JFIRST JFWLE Index of first x-grid point to be in- 

cluded in Cp output, JFWLE is the index 

of the first x-grid point on the airfoil 
in the extended fine grid 

JLAST JFWTE Index of last x-grid point to be include 

in Cp output, JFWTE is the index of 

the last x-grid point on the airfoil in 
the extended fine grid 

IBUF 0 Storage on logical unit 8 

0 - No storage 

1 - Unsteady results are being stored 

LIN 0 Selects formulas for Cp and local Mach 

number computations 

0 - Linearized expressions 

1 - Exact expressions 

IR r. 1 Program selected airfoil boundary 

1 condition based on IEQFLAG, Eq. (3) 

0 - No time-derivative terms 

(IEQFLAG = 0) 

1 - Include time-derivative terms 

(IEQFLAG = 1 or 2) 

DEL 6 .1 Airfoil maximum thickness to chord ratio 

XAR x 0. Location of airfoil pitch axis 









Table 2. Continued. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

ALFZRO 

a 

0 

0 . 

Mean angle of attack of airfoil (in 
degrees ) 

ALFONE 

“l 

1. 

Harmonic pitch amplitude of airfoil (in 
degrees ) 

XFLAP 

X f 

1. 

Control surface leading edge location 

XFR 

X h 

1. 

Control surface hinge location 

BETZRO 

*0 

0 . 

Mean control surface angle with respect 
to airfoil (in degrees) 

BETONE 

*1 

0 . 

Harmonic control surface amplitude (in 
degrees ) 

Z1 

Z 1 

0 . 

Plunge amplitude 

IMODFR 

r 2 

1 

Program selected wake condition based on 
IEQFLAG, Eq. (11) 

0 - No time-derivative terms (IEQFLAG=0) 

1 - Include time-derivative terms 

(IEQFLAG = 1 or 2) 

RSUB 


1.4 

Subsonic relaxation factor 

RSUP 


1.0 

Supersonic relaxation factor 

MITC 


500 

Maximum number of iterations 

MITM 


100 

on coarse, fine, and extended fine grids 

MITEF 


25 

for steady calculations, respectively 

TOLC 


.00001 

Convergence tolerance parameters in 

TOLM 


.000001 

steady calculations on course, fine 

TOLEF 


.000001 

and extended fine grids, respectively 

ICPRUN 


0 

Output request for unsteady pressures 

0 - No output 

1 - Output of unsteady pressures, if 

ISICO = 0 

LIND I F 


0 

Selects differential equation 

0 - TSD equation 

1 - TSD equation, nonlinear term omitted 
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Table 2. Concluded 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

ICPT 

r 3 

1 

Selects time derivative in linearized 
formulas for C p and local Mach number 

r.r- = 0 - Time derivative is set equal 
to zero 

r^ = 1 - Time derivative included in 
calculation of C p and local 
Mach number 











Table 3. Description of input parameters to NAMELIST MESH. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

JFMAX 


* 80 

Number of x-grid points in extended fine 
grid (Maximum value is 80) 

LFMAX 


61 

Number of z-grid points in extended fine 
grid (Maximum value is 61) 

XFINE 



Array of increasing x-coordi nates 
specifying the vertical lines in the 
extended fine grid (JFMAX elements) 

ZFINE 



Array of increasing z-coordi nates 
specifying the horizontal lines in the 
extended fine grid (LFMAX elements) 


*In case the default extended fine grid is used, it is sufficient to input 
the empty namelist $MESH $ 
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Table 4. Description of input parameters to NAMELIST AMPLI. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

N 

N 

1 

Number of terms in Eqs. (6) and (7) 

AB 

a »b 
n n 

1 

Array a n , b n in Eqs. (6) and (7) 
(N elements) 











Table 5. Description of input parameters to NAMELIST ORD. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

NINU 


★ 

Number of points on upper side of airfoil 

NINL 



Number of points on lower side of airfoil 

XINU 



Array of increasing x-coordi nates of 
points specifying upper side of airfoil 
(NINU elements) 

ZINU 



Array of corresponding z-coordi nates 
specifying upper side of airfoil (NINU 
elements ) 

VINU 



Array of corresponding slopes specifying 
upper side of airfoil (NINU elements) 

XINL 



Array of increasing x-coordi nates of 
points specifying lower side of airfoil 
(NINL elements) 

Z I ML 



Array of corresponding z-coordi nates 
specifying lower side of airfoil (NINL 
elements ) 

VINL 



Array of corresponding slopes specifying 
lower side of airfoil (NINL elements) 


* Default is the biconvex airfoil, 6 = .1 
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Table 6. Description of input parameters to NAMELIST COND. 


PARAMETER 

VARIABLE 

DEFAULT 

DESCRIPTION 

NIN 



Number of points at which chord de- 
formation g^(x) is specified, Eq. (5a) 

XIN 



Array of increasing x-coordinates of 
points where the chord deformation is 
specified (NIN elements) 

CONDI 

h 

0 

Array of corresponding g^values (NIN 
elements) 


C0ND2 


dg x 

dx 


0 


Array of corresponding dg^dx-values 
(NIN elements) 













Table 7. XTRAN2L extended fine default grid 


INDEX 

X 

z 

1 

-20.00000 

-25.00000 

2 

-16.30961 

-23.36111 

3 

-13.10054 

-21.77778 

4 

-10.33987 

-20.25000 

5 

-7.99453 

-18.77778 

6 

-6.03136 

-17.36111 

7 

-4.41705 

-16.00000 

8 

-3.11817 

-14.69444 

9 

-2.10110 

-13.44444 

10 

-1.33204 

-12.25000 

11 

-.77699 

- 11.11111 

12 

-.40170 

-10.02778 

13 

-.17160 

-9.00000 

14 

-.05175 

-8.02778 

15 

-.00667 

-7.11111 

16 

.00667 

-6.25000 

17 

.02000 

-5.44444 

18 

.04000 

-4.69444 

19 

.06000 

-4.00000 

20 

.08000 

-3.36111 

21 

.10000 

-2.77778 

22 

.12000 

-2.25000 

23 

.14000 

-1.77778 

24 

.16000 

-1.36111 

25 

.18000 

- 1.00000 

26 

.20000 

-.69444 

27 

.22000 

-.44444 

28 

.24000 

-.25000 

29 

.26000 

-.11111 

30 

.28000 

-.02778 

31 

.30000 

0.00000 

32 

.32000 

.02778 

33 

.34000 

.11111 

34 

.36000 

.25000 

35 

.38000 

.44444 

36 

.40000 

.69444 

37 

.42000 

1.00000 

38 

.44000 

1.36111 

39 

.46000 

1.77778 

40 

.48000 

2.25000 

41 

.50000 

2.77778 

42 

.52000 

3.36111 

43 

.54000 

4.00000 

44 

.56000 

4.69444 

45 

.58000 

5.44444 




Table 7. Concluded 


INDEX 

X 

z 

46 

.60000 

6.25000 

47 

.62000 

7.11111 

48 

.64000 

8.02778 

49 

.66000 

9.00000 

50 

.68000 

10.02778 

51 

.70000 

11.11111 

52 

.72000 

12.25000 

53 

.74000 

13.44444 

54 

.76000 

14.69444 

55 

.78000 

16.00000 

56 

80000 

17.36111 

57 

.82000 

18.77778 

58 

.84000 

20.25000 

59 

.86000 

21.77778 

60 

.88000 

23.36111 

61 

.90000 

25.00000 

62 

.92000 


63 

.94000 


64 

.96000 


65 

.98000 


66 

1.00000 


67 

1.02000 


68 

1.12274 


69 

1.35473 


70 

1.75323 


71 

2.35080 


72 

3.17695 


73 

4.25898 


74 

5.62249 


75 

7.29170 


76 

9.28970 


77 

11.63860 


78 

14.35968 


79 

17.47352 


80 

21.00000 





Appendix A - Results from Steady Calculations 
XTRAN2L Input 


NAC A 14 AO 1 0 A"fcB AlRFUll N4O.70O 4lFZ*0«l,0 BTtAOf.lTATt CONVANQCNCI 

SAEWOei t we 5P«5 » UtLTiO, 00*454369, UC4100.0, IfcMU 

J INPUT *N«0 t 780, XPK Au«0 , HA*JT«t02S NSHC*l28, ICMNT»120, [0UF*l, 
DCt»0«10» IARbO, 25, ALKZhO*| t OO# A17 QnEsq , 00 , UNO 
JHg3N A 
SAHPtt » 

SORQ NlMJ • 31# NJNL « Sir 
X JMj b 


•006667, 

.020000, 

.040000, 

,060000, 

,080000, 

•100000, 

• l 20000, 

,140000 » 

, 1 60000 1 

•160000, 

.200000, 

,220000, 

,240000, 

,260000, 

,200000, 

•JOOOQUr 

•320000, 

.34U000, 

,360000, 

,300000, 

,400000, 

,420000, 

,440000, 

,460000, 

,480000, 

,300000, 

,520000, 

,540000, 

,560000, 

,580000, 

,600000, 

,6200001 

, 6«P0nn, 

,660000, 

,600000, 

,700000, 

,720000, 

,740000, 

,760000, 

,780000, 

,800000, 

,960000, 

,820000, 

,440000, 

,840000, 

1,000000, 

,060000, 

,880000, 

,900000, 

,920000, 

,940000, 

ZZNU • 

•0)0)44, 

, 016833, 

.022381, 

,026732, 

.030912, 

,033687, 

•036532, 

,0391 !Sr 

,0«1413, 

,083883, 

,083258, 

.0*6879, 

,0*8300, 

,0*93J3, 

,050500, 

,051489# 

.052232, 

,03*732, 

,033113, 

,033200, 

.033113, 

,052*09, 

•US2SJ5, 

,051*21, 

,051033, 

.0*3333, 

,0**813, 

,0*7308, 

.0*6076, 

.0993**, 

,042865, 

.041114, 

»0392«0, 

.037367, 

,035365, 

.0133*3, 

,031231, 

,0**1*«, 

,026980, 

,024026, 

•0226*2, 

.005194, 

,020313, 

,002650, 

,018363, 

0,000000, 

.016182, 

,013*36, 

,011799, 

,009695, 

,007518, 

VI*U • 

,6663*6, 

,573066, 

,2*7303, 

.201161, 

.179388, 

,1303*7, 

,133839, 

,121231# 

,ioeoe*, 

,046760, 

,083*66, 

,076080, 

,06633*, 

,037031, 

,0983*7, 

,0*1197# 

.llJMO*, 

,022365, 

,009781, 

.,00050], 

«, 007*39, 

9,0137*0, 

8,0239*7, 

.,017319# 


•,056247, 

•,062229, 

.,068*31, 

.,078733, 

.,090909, 

.,0*33*0, 

.,089699# 

.,0*3739, 

•,09743b, 

.,100687, 

.,103*73, 

•.103608, 

.,106*02, 

107*39, 

.,1077*0# 

■>, 1 078 31 , 
-.1 itoit, 

•,108131, 

••130770, 

•,108983, 

-.133386, 

*,110036, 

-.109200, 

•, 1061*9, 

•>103139, 

..111830# 

XIM ■ 

,004467, 

,020000, 

,040000, 

,060000, 

,080000, 

.100000, 

,120000, 

,140000, 

.160000, 

,160000, 

,200000, 

,220000, 

.240000, 

•260000, 

,200000, 

,300000, 

.320000, 

,340000# 

,360000, 

,380000, 

,400000, 

,420000, 

,440000, 

,490000, 

.480000, 

,500000, 

,520000, 

,540000, 

,560000, 

.580000, 

,600000, 

,620000, 

.640000, 

,660000, 

,660000, 

, 700000, 

.720000, 

,740000, 

,760000, 

,780000, 

•eooooo, 

.960000, 

,620000, 

,480000, 

,840000, 

l.oooooo, 

,860000, 

,880000, 

,900000, 

,020000, 

,940000, 

ZINL ■ 

•,0089*l, 

•,016598, 

•,022527, 

-.026839, 

..030323, 

•,033748, 

•,036542, 

..03901*# 

-.0*1 .72, 

•,04325), 

.,0450)9, 

.,086684, 

• .089113, 

•.049J95, 

•,050481, 

.,031393# 

..052138, 

•,052706, 

•,05300), 

..033238, 

•,033286, 

•,053049, 

•,052654, 

;, 03*018# 

•,03120k, 

•,050192, 

•,049027, 

..0*7730, 

-.0*6317, 

•,044793, 

•,041)55. 

.,0*1*01# 

•,033538, 

•,037654, 

•,035699, 

.,033671, 

.. 011388, 

•,02*437, 

•,027263, 

.,0*3093# 

•,022831, 

•,00»»7S, 

•,020666, 

•,002450, 

•,018469, 

0,000000, 

.,0162*9, 

-.01397*, 

•,0)1722, 

•,009482, 

.,007*13# 

VI6L • 

..733823, 

•,364001, 

..**3107, 

•,1*91*6, 

•, 1*»*’0, 

•. l«7*21, 

•,134054, 

•,117213, 

•,103(33, 

•,092960, 

-,06*613, 

.,077611, 

•,06*016, 

9,03*1*7, 

•,04*8)2, 

•,04| J71, 

-.032330, 

•.02369J, 

.,013731, 

.,008037, 

.003**5, 

,01*397, 

,025)0), 

,036447, 

,038837, 

,054702, 

.061693, 

.0*7*7*, 

.073171, 

,079002, 

,064*07, 

.090100, 

,033308, 

,096371, 

,099283, 

,102*33, 

.106*33, 

.108179, 

,108874, 

,109377, 

,103838, 

.lt'228, 

,11046b, 

,122218, 

. 111*23, 
.123213, 

,112319, 

, 1 1*022, 

,112177, 

•112272, 

.114804, 


IEMO 

ICHND 0 



Appendix A - Concluded 
XTRAN2L Output 



PROGRAM ATRAN2L - VERSION 1.? 



illicit riNin difference program solving the complete 

UNSTEADY TRANSONIC SNAIL DISTURBANCE EOUATION 
INCLUDING AEROELASTIC RESPONSE 



UNSTEADY aerodynamics branch 



NASA LANGLEY RESEARCH CENTER 



HAMPTON* VA 23bfc 5 



CONTACT: DAVID A. SEIDEL CBO*> 865-A226 
JUhn T. BATINA (BOA) et&-A23c 




NACA MAO 1 0 AM|| AlRPQIL M«0,7IO AWFZMOrUO »TtAOY»STATI CONVERGENCE 

AtROlLAITXC SYSTEM 


IRESP 

5 

NO AX'POXl MOTION 

NELAST 

1 

2*00P , PLUNGE AND PITCH 


uc 

100,00000 

NONotHiNatoNAL velocity 

XHU 

10(00000 

HASS RATIO 


DElT 

,002u5 

TIME STEP , SECONDS 

IEOPLA0 

I 

Extran Equation 


anp l tO 

,01745 

AMPLITUDE op PULSE 

ROTHPLS 

10000(0 

PULSE EXPONENTIAL PACTOR 


tzerqst 

17,50000 

center op pulse 





INPUT PARAMETERS 






XX 

,20000 

REDUCED PNEQUCNCr BASED ON CHORD 

XM 

(TSOOO 

MEAN PREE-STREA" n *CH NO, 


GAM 

1,40000 

RATIO OP SPEC IP IC HEATS 





INEAO 

0 

initial conditions 

ICNVRG 

1 

GET S«S SOL, on OOARS, THEN PXNKft 

SNIPS 

MAX IT 

1024 

MAX, NU, TIME STEPS 

NSPC 

ia* 

NO, STEPS PER CTCLE 


ITRan 

0 

pitch, plunge, plap motion 

IS1C0 

0 

SINUSOIDAL VARIATION 


ICPPnT 

128 

ITERATION phint increment por CP 

I8UP 

1 

Output to lu o 


JPIN3T 

16 

PI»sr ChOMOmISE CRIOPT for cp PRINTING 

JL AS T 

6D 

LAST CHORORISE GHlDPT POR CP PRINTING 

LIN 

0 

LINEARIZED cp 

IR 

1 

TJHt-OHNIV, IN AINPCIL s,c. 


DEL 

,10000 

airfoil thickness ratio 

XAR 

,iaooO 

POSITION op PITCH Alls 


ALPZPO 

1,00000 

mean angle OP ATTACK 

alpone 

0,00000 

PITCHING amplituoi 


Z1 

0,00000 

PLUNGE AMPLITUOE 

IMODPK 

1 

T I N»*OEN I ¥ AT I It IN «AK( C0N0IT10N 


H SUB 

1,40000 

SURSUnIC RELAXATION PACTOR 

RSUP 

1,00000 

SuPtPSONJC PfLAIATlON PACTON 


NITC 

500 

maximum number op iteration* 

roue 

1 1 Q£*04 

tolmangi pakahitins in StCADT 


MJTM 

100 

ON COARSE, PINE AND EXT, PlNf GRID 

TOLM 

,10E«Qi 

CALCULATIONS UN COAttM, PINS 


MirtP 

25 

RESPECTIVELY 

fOLEP 

• 1 0C«03 

AND or, PIN! CHIU pcspectitilt 


ICPRUN 

0 

no REQUEST PQR OUTPUT OP UNSTEADY PRESSURES 

L I NO IP 

0 

NLHnTSP (9UAT10N («ITh NOn«LIn!AR 

TSPHJ. 

ICPT 

1 

TlME-OtR, COnTR J (J uTtS TO CP AND mACH, IP 1RU 




n \ number or (cqjsin harmonica 

MULTIPLIERS OP (CU)8IN HARH0HXC9 1,00000 


43 



Ntc* k*mo mi imoii ««o,rao alpino*i,o iTt«ur.»uTt coNvimNct 


XTCNAT1 ON NO, I 


ALP (DCS) • 

tint • 


1,0000000 

,0400171 


PLIINSX 
) (0(9) 


■ 0,0000000 

I 1,1119000 


inoia 


CPL 


*L 


OCPtCPl*CPU 


16 

,00667 

1,06346 

} ,36001 

,10730 

0,00000 


l 7 

,02000 

,2464] 

,33731 

,6104* 

,77817 


16 

,04000 

,53224 

,46393 

,76917 

• 71666 



,06000 

,25166 

.36042 

,10760 

,76172 


20 

,06000 

,20978 

,32207 

,62638 

,74061 


21 

,10000 

•17646 

•29272 

,63447 

,60312 


22 

, 1 2000 

•15361 

,26450 

,6*141 

,01302 


21 

,14000 

.12644 

,23042 

,84673 

, 0<2B] 


24 

,1*000 

•10166 

,21795 

,85440 

,62606 


25 

,16000 

,07916 

,19444 

,63661 

• *3221 


2<i 

,20000 

,05734 

•1*191 

,86229 

,63374 


27 

,22000 

,05373 

,13732 

,66636 

,04220 


16 

,24000 

,00901 

,12603 

,87102 

,65009 


it 

,26000 

••0142B 

,10121 

,87SOO 

•65619 


50 

,26000 

•,04030 

,07524 

,87407 

•66249 


51 

,30000 

•,06457 

,04007 

,16476 

,06613 


12 

,32000 

•,10406 

,01650 

,89233 

,67470 


53 

,34000 

•,13635 

•,01277 

,09931 

,68139 


12 

,36000 

•,16566 

••04J94 

,90266 

•98714 


53 

,56000 

•,16701 

-.07197 

,90331 

,84203 


36 

,40000 

•,21006 

••10120 

,90316 

,84669 


37 

,42000 

••24326 

•,14110 

,91061 

,90133 


36 

,64000 

••27841 

•,13923 

.91673 

,90494 


36 

,46000 

••30334 

••16150 

,91813 

,90379 


60 

,*6000 

••31680 

*,19747 

.91467 

,90398 


• 1 

,50000 

-,32716 

•,21034 

,91030 

,90107 


• 2 

,32000 

••33451 

•,22229 

,90713 

,04803 


«J 

,3*000 

•,35264 

•,23413 

,90463 

,69321 


44 

,36000 

•,3*451 

•,24630 

,90167 

,09274 


43 

,36000 

•, J7407 

•,23007 

,89002 

,69040 


46 

,60000 

••34223 

•,26716 

,69409 

,61703 


47 

,62000 

••38451 

•. 27271 

,64011 

,66244 


46 

,64000 

•,34366 

•,27673 

,86593 

,67734 


4« 

,66000 

-,40024 

•,26226 

,66142 

,67 J66 


30 

,66000 

•,40300 

• ,26093 

,67644 

,67062 


31 

,70000 

•,40303 

•,24333 

,67096 

,66679 


32 

,72000 

•,40204 

•,29301 

,86513 

,66146 


3 J 

,74000 

•,40099 

.,24194 

,63420 

,83334 


34 

,76000 

•,39426 

•,29059 

,63391 

,63027 


35 

,76000 

•,39607 

•,29044 

,0492} 

,64392 


36 

,60000 

•,39710 

•,29067 

,86492 

,64203 


37 

,82000 

•,39560 

•,24034 

,84069 

,63424 


36 

,84000 

•,39143 

•,26014 

,63532 

,83377 


36 

,06000 

•,31205 

• ,2»2U 

,62626 

,62715 


*0 

,88000 

*,36673 

*,27193 

,81967 

,62039 


61 

,40000 

••35671 

*,23661 

,81229 

,61144 


62 

,42000 

••3432? 

•,23960 

,10460 

,60123 


63 

,94000 

•,30636 

*,20306 

,78666 

,76342 


64 

,46000 

•,22266 

••14031 

,74346 

,75109 


63 

,46000 

*,06779 

•,03334 

,66164 

,69634 


66 

1,03000 

,36144 

,34040 

,31166 

,44617 

• , 

CQtPPXCIINTD 47 A NOLI 

OP ATTACK ■ 

1,0000000 

KlUN&C Q13P, • 

0,0000000 


L 

N6 M0M2NT 

(AIOUT I 6 

,23000 ) ■ 

•2330454C-O1 

NONHAl, PQ*CK • 

, 1 1643731*00 

A* J A 


,10107 

,00101 

1)110 

1071 

1109 

ilM 

1077 

1101 

1010 

2011 

2031 

2170 

1092 

l »*0 

1931 

1794 

till 

2370 

2202 

I960 

0044 

1217 

1014 

2140 

1012 

1442 

1721 

1473 

till 

1400 

1303 

1471 

140* 

1704 

1*07 

1041 

0000 

0903 

04*4 

074] 

044] 

0920 

10120 

,00044 

,00440 

,00410 

,10147 

,10124 

,04217 

,01423 

,04104 


,34410134*01 


44 




NACA 

66A010 A TIES 

AIRROII M*0 , 7 50 

AkFIKOFl.O .TEAOT..TATE CONVENGENCI 


iteration no, ■ 

1024 AIR. 

(OEQ) » 

1,0000000 

PLUNSI 

• 0,0000000 



TlMl 

• gi.ititiii 

t toes) 

■ 2660,0000000 


XNOtl 

X 

CPU 

CPU 

MU 

HI 

OCPOCPL-CPU 

it 

,00667 

,39640 

,97191 

,56009 

0,00000 

,61330 

17 

,02000 

• ,22203 

,36667 

,66676 

,56066 

,30090 > 

19 

,04000 

•,39426 

•12741 

,96147 

•71193 

,32167 


,06000 

•,61016 

,03017 

,96524 

.7641) 

,4411) 

20 

,06000 

•,44964 

••03271 

,98426 

,79663 

,41693 

21 

,10000 

•,47603 

•,07993 

,99497 

,02003 

,39600 

22 

,12000 

■,49301 

•,12362 

1,00211 

,64113 

,37019 

24 

,14000 

•,12607 

••16203 

1.01524 

.03923 

,36464 

2« 

,16000 

• ,514 1 0 

•• t’169 

1,02596 

,67303 

,36229 

23 

,16000 

• ,57309 

•,21605 

1,033)1 

,66407 

,33700 

2b 

,20000 

•,90709 

•,23519 

1,03674 

,69270 

,33190 

27 

,22000 

• ,60368 

••26091 

1,04520 

,90400 

,343J7 

29 

,24000 

■,62367 

-.28964 

1.05276 

,91691 

,33303 

29 

,26000 

• ,64170 

•,31540 

1,05960 

,92601 

,32630 

30 

,26000 

■,69414 

•,33537 

1 ,06449 

,93660 

,31917 

31 

,30000 

•,66657 

-, 35399 

1,06897 

,94433 

,31210 

32 

,32000 

•,68708 

•,37427 

1 ,07 663 

,95310 

,31282 

33 

,34000 

•,71842 

••394J8 

1,06823 

.96132 

,32404 

34 

,36000 

•,74912 

•,40942 

1 ,09946 

,96777 

,339/0 

31 

,36000 

-, 76851 

•,41792 

1,10653 

,97125 

,33019 

3ft 

,40000 

•,77«P7 

•.42273 

1,11027 

,97326 

, J56 1 4 

37 

,42000 

••79313 

•,42644 

1.11541 

,97479 

,36670 

38 

,44000 

• ,61940 

•,42490 

I.U961 

.97419 

,39442 

39 

,46000 

■,84691 

•,40970 

1,13444 

,96766 

,44664 

40 

,46000 

•,86263 

••36261 

1,14011 

,91660 

,47962 

•l 

,50000 

•,86757 

•,31319 

1.14184 

,944)7 

,31347 

42 

,52000 

•,65718 

•,32601 

1.13819 

,93260 

,33113 

43 

,14000 

• ,61066 

■,30019 

1,06296 

.92139 

,45007 

44 

,56000 

•,32049 

•,27742 

,93021 

,91147 

,04407 

41 

,58000 

•,19646 

25517 

.67514 

,90101 

•,01911 

46 

,60000 

-,20798 

-.23197 

,66040 

,89123 

•,02499 

47 

,62000 

•,2058? 

■,20147 

,67944 

,67926 

,00040 

4ft 

,64000 

•,19565 

17867 

,87477 

,66696 

,01696 

44 

,66000 

•, IBO07 

-,15522 

,66760 

,65603 

,02463 

10 

,66000 

•,16097 

-.13474 

,65674 

,64640 

,02623 

51 

,70000 

-, 1 J45 J 

•,11332 

,64866 

,03620 

,02619 

12 

,72000 

• ,11652 

•,06909 

,83774 

.62451 

,02743 

S3 

,74000 

.,09319 

p, 06383 

,62650 

,61214 

,02946 

5 4 

,76000 

•,07060 

•,0403b 

,61557 

,60047 

,03044 

51 

,78000 

• ,04964 

•,01910 

,60510 

,76975 

,03014 

16 

,80000 

•,02946 

,00097 

,79499 

,77950 

,03043 

17 

,82000 

• ,00985 

,02061 

,76504 

,76934 

,03045 

ia 

,64000 

,01077 

,04069 

,77445 

,73670 

,03012 

59 

,66000 

,03417 

,06)15 

,76203 

,74603 

,02050 

60 

,66000 

,06019 

,06723 

,74044 

,73301 

,02/04 

61 

,90000 

,08144 

• Mill 

,73097 

,72069 

,02967 

62 

,92000 

,09669 

,13470 

,72652 

,70736 

,03790 

63 

,94000 

,11607 

•16343 

,71677 

,69094 

,04146 

64 

,96000 

,16362 

•20616 

,09083 

,66170 

,04254 

61 

,98000 

,24943 

.27420 

,04162 

,62339 

,02477 

66 

1,00000 

,36469 

.3947J 

,14777 

,34037 

,01004 


ponce COEFFICIENTS AT ANGLE OF ATTACK • 1,0000000 FCUNOJ DIIF, * 0,0000000 

PITCHING FOMENT (AOOUT « ■ ,U000 ] • NQKHAt POKCC » »«I*l FUKCl ■ 


45 



Appendix B1 - Results from Unsteady Calculations for Harmonic Motion 

XTRAN2L Input 


N»c» 64*010 »H £5 *J»POIt Hao.710 *tPI»0«l,O AUPONEaO,* 14*0.071 
lUKOEl I’bSPao, uCa 1 30 a 0 1 <t«D 

JIsPlIT AKaO.11, KnaO.TBO, IkbAOaJ, M*KITH««o, XSPCaj60. ICPkkfaPOi IluFali 


0El«O,lO ( 

7 XA9«0,25# 

IlFZNOil 

.00* AU0N6iO ( 5Q, SEND 



i»t 3H S 
IAHPUI S 

lORQ NXNU 

a 51, njni 

, ■ 3 1 < 






XJNU ■ 

, 000007 , 

.020000, 

.340000, 

.060000, 

,060000, 

.100OQ0, 

,120000, 

,140000* 

.160000, 

.100000, 

.200000, 

,220000* 

.200000, 

.200000# 

,200000* 

,300000* 

.320000, 

.340000, 

.360000, 

,300000* 

.anoooo, 

,420000, 

,440000, 

,460000* 

,«§000n, 

.500000, 

.520000, 

,500000* 

.560000, 

.5*0000, 

,600000, 

,620000* 

, 0UOOOO , 

,660000, 

,690000, 

.700000, 

.720000, 

,740000, 

,760000, 

,700000* 

.800000, 

.960000* 

,020000, 

.900000, 

,940000, 

1,000000, 

,660000* 

.600000, 

.900000, 

,920000, 

,940000* 

ZJMM • 

.010166, 

,010499, 

.022191, 

,026732, 

,010412, 

.033687, 

,036532, 

,039113* 

.041413, 

.043444, 

.045259, 

.000074, 

.049100, 

,049535, 

,050504, 

,051403* 

.052232, 

.05279?, 

,053113, 

,033200, 

.013113, 

.0524(15, 

,052535, 

.051921, 

.051053, 

.049994, 

.049915, 

,007509, 

.046076, 

,044524, 

,042663, 

,041114* 

.01*230, 

,0371o7, 

.015395, 

,033303, 

.011251, 

.024124, 

,026960, 

,024020* 

.02267?, 

.0051*6, 

.020513, 

.002650, 

.018341, 

0 .000000, 

,010152, 

.013456, 

.011794, 

,009693, 

,007336* 

YlNU a 

.666JP», 

.37J0P6, 

.247501, 

.201161, 

.174596, 

.150597, 

,135030, 

•121230* 

.10808’, 

,045790, 

.095469, 

.076090, 

.066554, 

.057051, 

,040547, 

.041197, 

.0JJI04, 

.022195, 

.009743, 

-.000301, 

..00763’, 

..013740, 

••U2J947, 

-,037534* 

••049629, 

..056247, 

..062224, 

..069431, 

..0747 J3, 

..090404, 

-,005120, 

-.069094* 

..0’3M», 

-.097436, 

..100687, 

..103473, 

..105608, 

..1064o2, 

•,107439, 

-.107720* 

..107931, 

*.122939, 

.,109131, 

-.130770, 

. .108»91, 
-.111396, 

-.110036, 

-.104200, 

-.106129, 

..105359, 

-.111450* 

XINL • 

.006667, 

,020000, 

,0«0000, 

.060000, 

,080000, 

.100000, 

,120000, 

,140000* 

.160000, 

.100000, 

,200000* 

,220000, 

.240000, 

.260000, 

,200000, 

,300000* 

.320000, 

.340000, 

,360000* 

,300000, 

, 000000 , 

.420000, 

,440000, 

,400000* 

.430000, 

.500000, 

,320000* 

.500000, 

.560000, 

.560000, 

,600000, 

,020000* 

,640000, 

.660000, 

,600000* 

,700000, 

.720000, 

.740000, 

,760000, 

,760000* 

.000000, 

.960000, 

,820000, 

.900000, 

, 040000* 
} ,000000* 

.060000, 

.000000, 

.900000, 

,920000, 

,940000* 

HNL • 

.,009441 , 

.,016949, 

•,022527, 

•,0?b»39, 

••030529, 

••033740, 

•,030342* 

-, 039014* 

..041172, 

-. 043251 , 

-,0«5019, 

•,000600# 

• .naSU3* 

•,049395, 

•,050403, 

..03139J* 

..052139, 

.,092706, 

• .053001 * 

•.053256, 

-.053206, 

033049, 

-,052034, 

-.032036* 

-.051206, 

.,050142, 

• , 0«9027, 

•,007730, 

-.006317, 

•.004793, 

-,043153, 

•• 041401# 

. .039559, 

.,037654, 

•,03569a, 

•,033070, 

-.011500, 

•.029437, 

-.027205, 

•*025063* 

..022991, 

-.004975, 

-.020694, 

.,002459, 

• * 01 0609, 

o.oooooo* 

•,010229, 

-.013974, 

•.011722, 

-,009402, 

-.007215* 

VI ML « 

-.755623, 

.,364001 , 

•.263107, 

-.144126, 

-.109990, 

• a 1476?3, 

-.134034, 

-•117215* 

-.105643, 

-.092460, 

-.0*6615, 

.,077619, 

-.009010* 

-.059127, 

-, 049012, 

••041371, 

-.032*50, 

.,023643, 

-,013753* 

-.006057, 

.005205, 

•014597, 

,025101, 

,030447* 

.046437, 

,054702, 

•061095* 

.097974, 

.073371, 

•079002, 

,004907* 

,090160* 

.0*3900, 

,046172, 

,099265* 

.102431, 

.100255, 

.100179, 

,100074, 

,109377* 

.109836, 

•119226, 

.110406, 

.122214, 

, 1 1 1 025* 
,123213* 

.112514, 

. 112«22# 

•112177, 

,112272, 

,114604* 

SEND 
ICOnD s 
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Appendix B1 - Continued 


XTRAN2L Output 


naCA 64AQIQ knit AIRFOIL m*0 ,7 60 


AER0CLA8TXC SYSTEM 

XRE3P 

0 

HARMONIC OSCILLATION 

UC 

t 00 .00000 

n&nOImEnSJONAL velocity 

OELT 

.00116 

time step, seconds 

AMPtTD 

.017*5 

amplitude of pulse 

TZCPQ9T 

17,50000 

center of pulse 

INPUT PARAMETERS 
XK .19000 

REDUCED FREQUENCY BASED ON CHORD 

QAM 

1.40000 

RATIO OF SPECIFIC HEATS 

I READ 

5 

READ INITIAL CONDITIONS 

MAXIT 

1440 

MAX. NO. TIME STEPS 

ITRAN 

0 

PITCH, PLUNGE, FLAP MOTION 

ICPRnT 

90 

iteration print increment for cp 

JPIRST 

16 

first ChURdmISE GRlCiPT FOR CP PRINTING 

UN 

0 

LINEARIZED CP 

DEI 

.10000 

AlwFUJL THICKNESS RATIO 

AIPZRO 

1,00000 

-UN ANGLE UF ATTACK 

Z1 

o.ooono 

PLUNGE AMPLITUDE 

RSUB 

l .40000 

8U8S0NJC RELAXATION FACTOR 

MITC 

500 

MAXIMUM NUMgtfi OF ITERATIONS 

MI TH 

loo 

ON COARSE, FINE AND E * T , FfNC GRID 

Ml T£P 

2 5 

respectively 

ICPRUN 

0 

no request for output op unsteady pressures 

ICPT 

1 

TIME-DER, CONTRIBUTES TO CP ANO MACH, IF IR 

N 

l 

NUMBER of (CO) BIN HARMONICS 

"utTiftitRS or (cant** i.ooooo 


AlFZRQM 

•0 ALF0NE«0,5 

00,075 


NELAST 

l 

J-DOF 

, PLUNGE ANO 

PITCH 

AMU 

60,00000 

MASS RATIO 


IEGFLAS 

J 

Extran 

EQUATION 


mDThPLS 

10000,0 

PULSE 

EXPONENTIAL 

FACTOR 


AM 

,76000 

MEAN FREE-STREAM MACH NO. 


ICNVRQ 

1 

QET 8-6 SOL. UN C0AR8, THEN FINER 

srioj 

N3PC 

360 

no, STEPS PER CYCLE 


I8IC0 

0 

SINUSOIOAL VARIATION 


IBUF 

1 

Output to lu • 


JL AST 

66 

iter choro«ue I opt for cf frutino 

IR 

1 

TIME-DERIV, in AIRFOIL S.C. 


XAR 

,25000 

position of pitch axis 


alfune 

,50000 

PITCHING AMPLITUOE 


IMOOFR 

1 

Tihe-OERXVATIVC in make CONDITION 


R SUP 

1,00000 

SUPERSONIC RELAXATION FACTOR 


TQLC 

,10E*04 

tolerance parameters in steaot 


TOLM 

, l 0E-05 

CALCULATIONS on COARSE, FINE 


tolef 

.ioe-o3 

and ext, fine grid RESPECTIVELY 


L I NO IF 

0 

NLR-TSP EQUATION 1"XTn NDN-LINEAR 

TERMJ 


47 




N AC A 

64 AO 1 0 AMES 

AIRFOIL MR0.760 

ALPZROil.O 

aroNg.o.s k'0,071 


ITERATION NO, ■ 

l alf 

COCOS 9 

1,0067262 

PLUNGE 

• 

0,0000000 


TIME 

• 

10,2829357 

T (DEC) 

■ 

2661,0000000 


INDEX 

X 

CPU 

CPL 

MU 


"l 

DCPRCPL-CPU 

16 

,00667 

,35667 

,97144 

,56761 


0,00000 

,9l«57 

M 

.02000 

-.22173 

,36656 

,86669 


.56079 

,56629 

10 

,04000 

-,39406 

.12722 

,96144 


,71156 

,52128 

H 

,06000 

-.41042 

,03043 

,96624 


,76413 

,44064 

ao 

,04000 

44954 

-.03281 

,98430 


,79661 

,41671 

21 

,10000 

-, 47598 

-.06000 

,99499 


,62002 

,39367 

22 

.12000 

-.49381 

-, U363 

j, 00214 


,64106 

,37016 

2 J 

,14000 

52691 

-.16196 

1,01529 


,65918 

.3649] 

24 

, 1 b 000 

-, 55427 

-.19179 

1,02602 


,67296 

.36246 

25 

.18000 

-.57320 

-, 21590 

1,03337 


,66400 

.35730 

26 

.20000 

-.54729 

-.23499 

1,03681 


.69263 

,35230 

27 

.22000 

-.60414 

-.26025 

1.04527 


.90392 

.34389 

20 

.24000 

-.62396 

-.28953 

1,05263 


.91662 

,33446 

20 

.26000 

-, 64207 

-, 31503 

1,0596® 


,92792 

,32704 

30 

.26000 

-.45497 

-.33495 

1.06453 


.93651 

.32002 

31 

.30300 

-.66705 

-, 35350 

1 ,06905 


.94443 

,31356 

32 

.32000 

•.6876J 

-.37372 

1,07672 


.95100 

,31191 

33 

.34000 

-.71902 

-.39377 

1,08633 


,96141 

,32525 

34 

.36000 

-. 74979 

-.40875 

1,09958 


.96766 

.14104 

33 

,38000 

-.76924 

-.41719 

1.10662 


.97117 

,35205 

36 

.40000 

-,779fcfc 

-.42194 

1. U037 


.97315 

,35772 

37 

.42300 

79399 

-, 42558 

1,11551 


.97467 

.36640 

3» 

,44000 

-.82031 

• e 42«06 

1,12491 


,97407 

.39625 

30 

,46000 

-.84751 

-.40873 

1.13455 


,96776 

,43660 

aO 

,48000 

-.66367 

-.38177 

1,14022 


.95655 

,46191 

«1 

.50000 

*, 86867 

35249 

1.14195 


,94423 

,51618 

42 

,52000 

-.85635 

-.32486 

1.13631 


,9324b 

,33347 

aj 

.54300 

-, 65189 

-.29915 

1.06311 


.92145 

,35254 

aa 

.56000 

-.32179 

-.27612 

,93036 


.91132 

,04567 

45 

.58000 

-, 19782 

-,25420 

,67530 


.90165 

-,05936 

46 

.60030 

-.20941 

-.23054 

,66056 


,89109 

-.02113 

47 

.62000 

-.20736 

-.20398 

,87961 


,87909 

,00339 

40 

.64000 

-, 19721 

-.17711 

,87494 


,66678 

,02010 

49 

,66000 

-,16173 

-, 15360 

,86778 


.85587 

,02610 

50 

.68000 

-, 16267 

-.13305 

,85892 


.64622 

,02962 

51 

.73000 

-.14127 

-.11157 

,64665 


.63602 

,02970 

52 

.72000 

*. 11834 

-.06727 

.63793 


,82431 

.03107 

53 

.74000 

-.09406 

-,0b|95 

,82669 


,81194 

,01114 

54 

.76000 

-.07275 

-.03640 

,81577 


,8002b 

,03435 

55 

.78000 

-,05166 

-.01706 

,80531 


,76954 

,03456 

56 

.80000 

-,03155 

.00306 

,79521 


.7792B 

,03460 

57 

.62000 

-.01200 

.02276 

,78526 


,76911 

,03476 

50 

,84000 

,00855 

.04511 

,77467 


,7584? 

,03456 

50 

,86300 

.03228 

.06544 

,76226 


,74661 

.03315 

60 

.66000 

.05764 

.06956 

,74866 


,73357 

,03174 

61 

.93000 

,07902 

.11353 

,73722 


,72040 

,03451 

62 

.92000 

.09443 

.13727 

,72877 


,70710 

,04267 

63 

,94000 

.11552 

,16596 

,71703 


,69067 

,05047 

64 

,96000 

,16101 

.20876 

,69111 


,66541 

,04777 

65 

.94000 

,24274 

,27669 

,64192 


,62308 

,03414 

66 

1.00000 

.38359 

.39604 

,54692 


,54106 

.01245 

FORCE COEFFICIENTS 

AT ANGIE OP 

ATTACK 4 

1 ,0087261 

PLUNGE OZSP, 

• 

0,0000000 


AIRFOIL 

PITCHING MOMENT (ABOUT X ■ ,25000 ) 1 

. l 0 36569E-01 

NORMAL FORCE 

• 

,234il69E»0Q 

AXIAL FORCE ■ 


,»]«0)t»C0i 
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N AC A 

AAA010 AMES 

Al»POil H*0 ,7 60 

ALP2R0- 1 ,0 ALPONftO.S 4*0,075 


ITERATION NO, « 

U40 ALP 

(DEC) • 

1,0000000 

PlUN&e ■ 

0,0000000 


TINE 

1 T1.J9S22JT 

7 (DEO) • 

4320,0000000 


inqe* 

X 

CPU 

CPI 

MU 

"U 

OCPiCPl»CPU 

i* 

,00667 

,40822 

,92937 

,33056 

0,00000 

,52114 

17 

,02000 

••16362 

,32377 

,66007 

•59123 

,48739 

16 

,04000 

•,34306 

,09094 

,93987 

,73235 

,43402 

19 

.06000 

•,36648 

•,00023 

,95055 

,76074 

,36825 

20 

,06000 

•,40829 

■,06062 

,96712 

,61122 

, J47 66 

21 

,10000 

•,43818 

•,10590 

,97933 

,63332 

,33238 

22 

,12000 

-, 45974 

..14796 

,96807 

,65341 

,31179 

2) 

,14000 

•,49128 

•,18506 

1,00072 

,67074 

,30620 

24 

,16000 

•,52130 

•,21369 

1,01262 

,68366 

,30761 

29 

,16000 

• •34240 

•,23666 

t ,02087 

,69427 

,30574 

20 

,20000 

•,5540! 

•,25459 

1,02537 

,90233 

,29942 

27 

,22000 

• ,56993 

27921 

1,03152 

,91324 

,2907| 

2d 

,24000 

•,39210 

•,30824 

1,04005 

,92595 

,28386 

29 

,26000 

••61330 

■,33348 

1.04814 

,93685 

.27983 

SO 

.26000 

••62733 

•,35291 

1 ,05344 

.94518 

,27442 

31 

,30000 

-.63418 

..37103 

1,05790 

,95287 

,26816 

32 

,32000 

••65974 

-.39111 

1,06561 

,96113 

,26863 

S3 

.34000 

•,69167 

‘ •,«U23 

1,07756 

.96973 

,28065 

34 

,36000 

• ,72333 

•,42596 

1,08917 

,97584 

,29737 

S3 

.39000 

-.74276 

•,43355 

1 ,09627 

,97900 

,10924 

SO 

.40000 

••7527? 

•.43718 

1,09986 

,98053 

.31534 

J 7 

.42000 

•,76667 

•.43962 

1,10466 

,98166 

,32685 

36 

,44000 

•,79298 

•,43696 

1.11434 

,96053 

,35602 

19 

,46000 

-.61966 

•.41945 

1,12166 

.97318 

,40023 

40 

,48000 

-.63306 

•,39o30 

1,12858 

,96134 

,44278 

41 

,50000 

•.82691 

•,35952 

1.12704 

, 9u64 5 

.46940 

• 2 

,52000 

•,70356 

•,33090 

1 ,08144 

,93630 

,37266 

43 

.54000 

•.41870 

•,30463 

,96943 

.92502 

,11*06 

44 

,56300 

-.25661 

•,28065 

,69996 

,91466 

•,02403 

43 

,59000 

•,26264 

•,25649 

,90266 

,90485 

,00415 

46 

,60000 

•.25718 

•,23443 

,90029 

,69414 

,02276 

47 

.62000 

•.24129 

-.20750 

,69413 

,66201 

,03579 

46 

.64000 

•.22479 

••16032 

,66562 

,66959 

,04447 

49 

.66000 

••20335 

• , 1 Sb5 9 

,87609 

,65660 

,04675 

SO 

.66000 

••18001 

•,13589 

,66536 

,6u869 

,04412 

SI 

,70000 . 

15542 

• ,1U27 

.85391 

,83664 

,04115 

32 

.72000 

-.1300* 

•,06965 

,84193 

,82690 

,04022 

S3 

, 7 “ 000 

10491 

•,06442 

,62968 

,81449 

,04030 

34 

,76000 

•,06110 

-,04062 

,81630 

,80280 

,04028 

55 

,78000 

•,05679 

•,01947 

,60730 

,79206 

, 039 J2 

S6 

,60000 

•,03766 

,00066 

,79674 

,78164 

,0383 i 

57 

.62000 

•,01724 

,02034 

,78641 

,77169 

,03738 

S6 

,64000 

,00408 

,04064 

,77547 

,76106 

,03637 

59 

,66000 

,02851 

,06291 

,7627* 

,74927 

,03441 

60 

,68000 

,05469 

,08696 

,74886 

,73629 

,03229 

61 

,90000 

,07646 

,11063 

,73712 

,72320 

,03437 

62 

,92000 

,09240 

,13443 

,72640 

.71001 

,04201 

63 

,44000 

,11«07 

,16299 

,71639 

,69370 

,04892 

64 

,96000 

, 1 601)6 

,20561 

,69017 

,66863 

,04533 

65 

,96000 

,24226 

,27151 

,64066 

,62661 

.03123 

46 

l.ooooo 

,36270 

,39366 

,54560 

,54422 

,01116 

FORCE COEFFICIENTS AT iNELt OF 

ATTACK ■ 

1,8000000 

PIUNCC OSIP, • 

0,0000000 


AIRFOIL 

pITCmJnC moment 

(ABOUT 1 • .29000 ) • 

,10497b Jt*01 

NORMAL PORCE • 

.2013339E400 

4IIAL 70RCI ■ •.3739490E-02 
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N6C* 14*010 4MES aIRFOJL "40.T80 *IFI»0«1,0 


*lFQNt»O.J K10.075 


LAST OOCILIATOHY CYCLE OF CF TIfe FlSTOFY 


0105? 

,01063 

,01070 

,01076 

,01083 

,01069 

,01096 

,01103 

,01109 

•01115 

01122 

,01129 

.01135 

.01141 

,01146 

,01151 

,01160 

,01166 

,01173 

.01179 

ones 

,01191 

,0119? 

.01203 

,01209 

,01215 

,01221 

•01227 

,01(33 

•01239 

01244 

,01250 

.01256 

,01261 

,01267 

,01273 

,01276 

,01164 

,01269 

,01294 

0129* 

.01305 

,01310 

,01315 

.01320 

.01325 

,01329 

,01334 

,01339 

,01343 

oisua 

,01352 

,01357 

.01161 

.01365 

.01369 

,01373 

,01377 

,01361 

,01365 

01369 

,01391 

,01396 

.01199 

.01402 

,01405 

,01409 

,01412 

,01414 

• 01417 

0J42Q 

.01422 

•Ol«25 

.01427 

.01429 

.01432 

, 0 1 4 J4 

.01435 

,01437 

.01439 

01««0 

,01492 

,01443 

,01444 

.01445 

,01446 

.01447 

,01446 

,01448 

.01449 

01049 

,01449 

,01449 

.01449 

.01449 

.01448 

,01446 

.01447 

.01446 

,01446 

uiooo 

,01441 

,0144| 

,01440 

,01439 

.01437 

,014)5 

,01432 

,01430 

,01428 

0 1 «25 

,0I«22 

,01«19 

.01416 

.0H13 

,01409 

,01406 

.01402 

.01346 

.01394 

01390 

•01385 

,01361 

.01376 

.01371 

,01366 

,01361 

.01155 

,01350 

.01344 

01336 

,01332 

.01326 

,01119 

.01311 

.013U6 

,01299 

,01292 

,01265 

.01276 

012T0 

,01262 

,01255 

.01247 

,01219 

.01231 

,01222 

,01ll« 

,01205 

,01197 

01186 

,01179 

,01170 

.01161 

.01152 

,01143 

,01133 

,01124 

,01114 

,01104 

01095 

,01065 

,01075 

,01065 

,01055 

,01045 

.01035 

,01025 

,01015 

•01004 

00994 

,00964 

,00974 

,00961 

,00953 

,0094 J 

,00932 

,00922 

,00912 

•00901 

00192 

,00101 

,00671 

,00461 

,00651 

,00641 

,00631 

,00821 

,00611 

,00601 

00792 

,00762 

,00772 

,00761 

,00753 

,00744 

•007J3 

,00726 

,00717 

,00708 

00999 

•00690 

,00662 

,006? 4 

•00665 

,00657 

.00649 

,00642 

,00634 

•00627 

006 1 9 

,00612 
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Appendix B1 - Concluded 
Post-Processed Pressure Harmonics 


PRE88URE COEFFICIENT HARMQNlCS 
UPPER PRESSURE COEFFICIENT 


X 

■ *0047 


,0200 


,0400 


,0600 


,0600 


,1000 


• 1200 


,1«00 



real 

IHAC 

REAL 

IHAG 

real 

IHAG 

HEAL 

IHAG 

>(»L 

IHAG 

REAL 

1HA6 

REAL 

IHAG 

REAL 

l H A 8 

0 

.3363 
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3.9173 
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Appendix B2 - Results from Unsteady Calculations for Pulse Motion 

XTRAN2L Input 
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Appendix B2 - Continued 
XTRAN2L Output 
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Appendix B2 - Concluded 
Post-Processed Generalized Aerodynamic Forces 
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Appendix B3 - Results from Unsteady Calculations for Aeroelastic Motion 


XTRAN2L Input 


NACA 64A010 AMES 4 JMFUJL M*U,70O A^FZNClal ,0 U/C*3J0#Q AtKQtL A 8 T IC 
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Appendix B3 - Concluded 
XTRAN2L Output 


NACA AAA010 AMES AlMFOJL HSO,780 


AtrirtO « l ,0 


U/COJO.O 


AtROELASTXC 


AtRQEU A 3 TIC ST STEM 


1*EBP i 
UC 330*00000 
OCLT .00050 
AHPLTO ,01745 
TZEflOST 17,30000 
MALPHA 1,86500 
OHEGAA 100,00000 


aeroelastic transient 

NONOI^tNSIUNAL VELUC1TY 

TINE SUP* SECONOS 

AriPtPUOE OF PULSE 

CENTER OP PuLSE 

AIRFOIL RADIUS OP GYRATION 

UNCUUPLEO PITCH naT P*£U (RAO/SEC) 


NELAIT I 
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ItOPLAG > 

■ OThPLS 10000,0 
A ALPHA 1,80000 
UnEGAH 100,00000 
Ah *2,00000 


2. OOP t PLUNGE ANO PITCH 
HASS RATIO 
fcXTKAN EQUATION 

PULSE EXPONENTIAL factor 

UlST FROM ILAST AXIS TO AI KF HAI8 CENliN 

UNCOUPLEO PLUNGE NAT FME8 IRAO/SECJ 

elastic avis locatiun 


INPUT PARAHETERS 


XK 
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RATIO OF SPECIFIC HEATS 
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XM 
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ISICO 
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X A R 
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u 

1 

66 

1 
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l 
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0 
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SINUSOIDAL VARIATION 
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PITCHING AMPLITUDE 
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THE A MATRIX ISI 


4 « 

0 . 0 , 

o. o. o 

•,100006*05 0* o 

0. -,86956t*04 o 


THE B MATRIX 191 

4 2 

0 . 0 . 

0 , 0 , 

,100006*01 0. 

0. ,100006*01 


THE C MATRIX 181 

4 4 

,100006*01 0, 0 

0* • 1 000QE+01 o 

0 . 0 , 

0 # 0 , • 


f 1 00006*0 1 0. 

• ,10000E*01 

. 0 , 

• 0 , 


• 0 . 

• o • 

,100006*0} •,900C0E*00 
,900006*00 , 6b95bt*00 


matrix entering Invrhx i a a 

4 4 

,100006*01 0, 0. 0, 

0. • 1 00006*01 0. 0, 

0. 0, ,100006*01 .,900006*00 

0, 0, •,900006*00 ,869966*00 


INVERSE OF MATRIX 181 
4 4 

,100006*01 0, 0, 0, 

0, ,100006*01 0, 0, 

0, 0, ,146016*02 ,151126*02 

0, 0, ,151126*02 ,1679(6*02 


THE REOUCEO SYSTEM IS • 


THE A matrix 181 

4 4 


0 , 0 , 

0 , 0 , 

•,146016*06 -,131416*06 
• ,131 126*06 •,1460U*06 


,100006*01 0, 

0, ,100006*01 

0 . 0 . 

0 , 0 , 


THE M MATRIX ISI 


4 


2 


0 , 

0 , 

,146016*02 

,131126*02 


0 . 

0 . 

,151121*02 

,lb791t*02 


THE TRANSITION matrix 7 TERMS 19 I 
4 4 

,961066*00 •,16326t«01 ,«9697f03 •,27|76C-0l 

-, 167736.01 ,961666*00 •,31)666.09 ,496976.03 

• ,721466*02 •,649066*02 ,961646*00 •,U3266«01 

• ,7464)1*02 -,721«96*02 • 1 1 8 7731*0 1 ,9Sj66tt00 



THE Integral OF The Transition NATNII isi 


4 4 

,446476.01 -,27276t-03 ,124626-06 -,la|37E*09 
-.IIlbSE.OS ,496976-01 -.192386-09 , 12A62E-Q6 
.,l»|4«E.OT ., 161266-01 .44*471-01 -,272766-05 
-,167T5E-0l -,181446-01 -, 113686-03 ,496976-01 


THE FINAL HATNICE9 • 


THE PHI MATRIX 191 

4 * 

,48l86E*00 -, 163266-01 # 49b97C-03 -,272766*05 

-,i8775t-oi ,9Siabt*oo -.nibet-os , 496476. ui 

721481*02 -,6*9066*02 ,981866*00 -,161266-01 

.,746416*02 -,72I48t*02 -.187756-01 ,98llbt*00 


THE TH6TA*B-PHIHE MATRIX 181 

4 2 

,161446.05 ,187736-05 

,167736-05 ,208666-03 

,721466-02 ,746416-02 

,746416-02 .824716-02 
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(a) Pulse input. 


(b) Force output. 


(c) Unsteady force response 


Fig. 1 Pulse transfer-function analysis for unsteady aerodynamic forces; 







fa rf i eld boundary conditions. 
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(a) Normal force coefficient and pitching moment coefficient time 
histories. 

Fig. 5 Results from steady calculations (no airfoil motion, IRESP = 5) for 
the NACA 64A010A airfoil at M„ = 0.78 and oq = 1.0°; 
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a, deg 


(a) Four cycles of simple harmonic pitching motion about the 
quarter-chord. 


Fig. 6 Results from unsteady calculations due to harmonic pitching motion 
(IRESP = 0) for the NACA 64A010A airfoil at M„ = 0.78, a 0 = 1.0°, 
k = 0.075, and a x = 0.5°; 
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(b) Resulting normal force coefficie 
and reference sine wave. 


Fig. 6 



Reference 
sine wave 


pitching moment coefficient 




1 .2 i 


Upper 

Lower 



(c) Mean pressure distributions. 


Fig. 6 Continued. 
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(e) Real and imaginary components of the lower surface pressure 
coeff icent. 


Fig. 6 Continued. 
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(f) Real and imaginary components of the lifting pressure 
coefficient. 


Fig. 6 Concluded. 




(a) Pitch pulse input. 

Fig. 7 Results from unsteady calculations due to pitch pulse (IRESP 
the NACA 64A010A airfoil at M*, = 0.78 and oq = 1.0°; 


, deg 


3) for 


79 





Real 

Imaginary 



(a) Pitch and plunge aeroelastic transients. 


8 Results from unsteady calculations for aeroelastic motion (IRESP = 1) 
for the NACA 64A010A airfoil at M* = 0.78, oq = 1.0°, and 
U/c = 330.0; 




(b) Normal force coefficient and pitching moment coefficient 
aeroelastic transients. 


Fig. 8 Concluded. 
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